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_ otating motion is a design criterion of rotating 
machinery, and vibration, as a side effect, diverts 
„energy from this desired condition. There are 





_ | Agnes Muszynska, PhD three major categories of vibrations in mechanical sys- 

_ | Senior Research Scientist and tems: free, forced, and self-excited. These categories are 

| Research Manager based on fundamental mathematical models of these 

| Bently Rotor Dynamics Research vibrations, adequate to the observed physical phenome- 

E Corporation na. The main characteristics of these three vibration cat- 

| agnes@brdre.com egories are shown in Figure 1 and described in Table 1. 
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Figure 1. Three categories of vibrations: free, forced, and self-excited, and their characteristic responses, 
@, is a system natural frequency. 
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The major distinction between free, forced, and self- 
excited vibrations is the way in which the energy, nec- 
essary to sustain vibrations, is provided to the 
mechanical system. 


Free vibrations are excited by an 
impulse or step force applied to the sys- 
tem and occur with one or several system 
natural frequencies. In rotating machines, these impuls- 
es may be due to an electrical short circuit in a motor or 
generator, a rub, a surge, a slug of water, a broken blade, 
etc. If the resulting response is not so high as to cause 
internal damage, and if the system is stable (has no 
internal sources of energy equal to or greater than its 
losses), then the impulse-excited free vibration will 
decay as the initially input energy dissipates, at a rate 
which depends on the amount of system damping. The 
adequate models of the mechanical system free vibra- 
tions are usually presented as linear, autonomous (no 
external, time-dependent excitation) ordinary differen- 
tial equations. 





An unstable mechanical system (containing an inter- 
nal source of energy greater than its losses) will respond 
with increasing amplitude vibration. Linear models 
(see Appendix 1) of such a system predict only the insta- 
bility threshold, and an unrealistic, infinite growth of 
vibration amplitudes in the post-instability-threshold 
domain. Practically, therefore, the linear models 
become inadequate when the system becomes unstable. 


The forced vibrations are excited by 
continuously applied, time-dependent 
external forces. Most often these are peri- 
odic forces. The most common periodic 
force in rotating machines is unbalance. 
Other, nonsynchronous, forces could be caused by vane 
passing, element passing in rolling element bearings, 
gear tooth mesh, etc. The mechanical system response 









follows the excitation, with some phase lag. The fre- 
quency of forced vibrations is the same as the excitation 
frequency (provided that the system is linear). A reso- 
nance occurs when an excitation frequency coincides 
with one of the natural frequencies of the system. Fora 
mechanical resonance, the forced vibration amplitude 
peaks, and the phase of the exciting force and the lag- 
ging vibration response differ by 90 degrees. 


The adequate models describing forced vibrations of 
mechanical systems are linear, non-autonomous (with 
time-dependent excitation), ordinary differential equa- 
tions. They are similar to free vibration models, except 
for the addition of the exciting force (the non- 
autonomous term). 


The self-excited vibrations belong to the 

Ge third category. The energy supply to sus- 
"So % tain the periodic vibrations is constant in 
this case, and comes from a source which 
may be internal or external to the system. 


The system provides the energy dosage from the con- 
stant supply through an internal energy converter, as 
part of a feedback loop, with a frequency characteristic 
to it; that is, one of the natural frequencies of the system. 
Self-excited vibrations of machinery rotors can take the 
form of fluid-induced (oil, steam, or gas) or rotor inter- 
nal (material) hysteresis-induced whirl or whip. The 
energy to sustain these vibrations is often drawn from 
the rotational motion of the rotor. 


Adequate models of self-excited vibrations, because 
they must include nonlinear terms, are nonlinear, ordi- 
nary differential equations. Such nonlinear models very 
often represent the extension of the linear autonomous 
models for the systems which are potentially unstable, 
due to the existence of an internal source of energy. The 
addition of the nonlinear terms to these models elimi- 
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Table 1. Characteristics and features of the three types of vibration 
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nates the incorrect prediction that, just after the instabil- 
ity threshold, vibration amplitude grows infinitely. As 
vibration amplitude increases, the nonlinear term 
becomes dominant and slows down the amplitude 
growth. The vibration amplitude continues to grow 
until it reaches the limit cycle of the self-excited vibra- 
tions. The self-excited vibrations are then sustained by 
the constant energy source, The limit cycle is a reflec- 
tion of the balance between the linear and nonlinear ele- 
ments in the system. 


For fluid whirl or whip, the constant source of ener- 
gy is shaft rotation, and the energy converter is fluid 
dragged into circumferential motion by friction. 


Free and forced vibrations: linear model 


The mathematical model of the one forward lateral 
mode isotropic rotor is: 


Miz + D+ Kz+ D(z—jAgQz)= Fata) 


“= d/dt (1) 


z=x+jy, j=v-l, 


where x(t), yf) are rotor lateral displacements lumped 
into one lateral coordinate 2(r); M, D., K are rotor modal 
mass, external damping, and lateral stiffness, respec- 
tively; D, Aare fluid environment damping and fluid cir- 
cumferential average velocity ratio, respectively; 2 is 
rotor rotative speed, which represents an internal source 
of energy of the rotor system; F, @, 6 are amplitude, 
frequency and angular orientation of the nonsynchro- 
nously rotating force applied to the rotor; and E is time. 


The rotor model (1) is a non-autonomous, linear dif- 
ferential equation with the complex lateral coordinate 


zit). If the external exciting force does not exist, (F=0), 
then the model becomes autonomous. 


The autonomous (F=0) equation (1) predicts the 
threshold of instability [1] 


l DTK 
Q.=—)14+—£|_/— (2) 
e 


where 42, is the rotative speed at which the free vibra- 
tions become harmonic {rotor circular orbit): 


A) = ce © 


with the natural frequency VK/M and a constant 
amplitude C, as the damping force becomes nullified by 
the action of the tangential force j/DAQz. For rotative 
speeds greater than 42, the linear model predicts an 
exponential growth of free vibration amplitudes: 


z(t) = Celtis } (4) 
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where œ, is the “damped” natural frequency of the sys- 
tem, and y (Appendix 2) is a positive number. 


In this situation, the linear model (1) is inadequate to 
describe the behavior of the rotor; it needs to be com- 
plemented by the nonlinear terms which will be dis- 
cussed in the next section. 


5o far only the autonomous part of the equation (1) 
was considered (F=0). Now the concentration will be on 
forced (or “excited”) vibrations of the system, in 
response to the periodic, nonsynchronous exciting force 
on the right side of Eq. (1). 

The rotor forced response will be periodic with the 
same frequency as the exciting force: 


z(t) = Aea (5) 


where A and œ are the response amplitude and phase, 
respectively. 


F 
(x -Mo 4 (D+ D)o- paf ° (6) 





(D+ D, jo- DAQ 
a = 6 — arctan == 
K = Mor 


The peak amplitude, referred to as a resonance, will 
occur at the frequencies @ at which either the Direct 
Dynamic Stiffness [1] K - Mæ? is null (in the case of low 
damping) or the Quadrature Dynamic Stiffness 
(D + DJ@- DAQ is null (in the case of high damping). 
In the first case, the resonance is referred to as a 


“mechanical resonance,” and in the second case as a 
“fluid-induced resonance.” In the first case, the reso- 


nance frequency is © = v K/M | in the second case @ = 
AQ2/(1+D_/D) , both of which represent rotor natural fre- 
quencies. 


Self-excited vibration: nonlinear model 


The nonlinear model of the one forward lateral 
mode isotropic rotor results from Eq. (1) by the addition 
of anonlinear term. For simplicity, assume that only the 
rotor radial stiffness is nonlinear, and take into account 
only the first term of the power series of the isotropic 
nonlinear stiffness function [2]. The rotor nonlinear 
model is, therefore, as follows: 

Mi+ D ¿+ Kz+ D(z jaQz)+ K,|z/'2=0 (7) 


izl i x + yi 
where K, is the coefficient of the nonlinear stiffness term. 


The external exciting force is omitted in Eq. (7) for clari- 
ty, resulting in an autonomous system. 
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It will be shown below that Eq. (7) has an exact solu- 
tion describing the limit cycle self-excited vibrations 
(Figure 2): 


a(i) = Be’ i 


where œ, B are the frequency and amplitude of the self- 


excited limit cycle, which can be calculated when Eq. (8) 
is substituted into Eq. (7) 


|-Mw* + jD,w, + K+ jD(w, -AQ)+ KB |Be™ =0 (9) 


To satisfy this equation, the expression inside the 
brackets must equal zero. ‘Splitting it into real and 
imaginary parts, the unknown frequency @, and ampli- 
tude B can be calculated: 

a (10) 


The transient process, which starts at the instability 
threshold and ends at the limit cycle, can be calculated 
when the solution of Eq. (7) is assumed as follows [3]: 


z=y(rje™ (11) 


an MEA 
( 


ü, = —— a, B= Miaka -= 
l+ D/D | K, | (1+ D/D) 


where wr) represents a variable response amplitude, 
and is the frequency at the instability threshold. For 


the model (7) @,=VA/M, 


Oil whirl limit cyete of 
sell-exciled vibralbons 


Small 1% vibration just 
before the Instability Threshold 





Figure 2. Self-excited vibration transient process: increasing 
rotor vibration amplitude leading to the limit cycle of oil 
whirl, Each white dot in this oscilloscope picture corre- 
sponds to one rotor rotation. 
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Appendix 1 


Linear terms in an ordinary differential equation con- 
tain only terms with the first powers of acceleration (z) 
velocity (z), and displacement (z) with some constant (or 
time t dependent) coefficients. If z, z, Z appear in the 
equation as Sells, higher powers, or products (for 
example, sin(z), z, or zz), the differential equation 
becomes nonlinear. 


Appendix 2 
The solution of the eigenvalue problem for Eq. (1) 
provides two eigenvalues s, , [1]: 
(Pe cer Pee 
| ) 
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where @ is the “damped” natural frequency, 
| 42 
oe | 
= E + ; 
ii ala \ í M, 
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M 2M 











and the real part of one eigenvalue (A1) is always neg- 
ative. 


The instability threshold occurs when the real part of 
the other eigenvalue (A1) becomes zero. The instability 
condition (2) results from the solution of the second real 
part set equal to zero: 


pp ee (282). 
M 


ae i er aE JE + 
2M 2 Y : 
At the instability a the natural frequency 


becomes “undamped” and is equal to vE/M, For the 
post-instability-threshold rotative speeds, the real part 
of one eigenvalue becomes positive (see Eq. (4)): 
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as the term (DAQ/M) becomes dominant. 
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